Length requirements for numerical-relativity waveforms 
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One way to produce complete inspiral-merger-ringdown gravitational waveforms from black-hole-binary sys- 
tems is to connect post-Newtonian (PN) and numerical-relativity (NR) results to create "hybrid" waveforms. 
Hybrid waveforms are central to the construction of some phenomenological models for GW search templates, 
and for tests of GW search pipelines. The dominant error source in hybrid waveforms arises from the PN contri- 
bution, and can be reduced by increasing the number of NR GW cycles that are included in the hybrid. Hybrid 
waveforms are considered sufficiently accurate for GW detection if their mismatch error is below 3% (i.e., a 
fitting factor above 0.97). We address the question of the length requirements of NR waveforms such that the 
final hybrid waveforms meet this requirement, considering nonspinning binaries with q = M2/-M1 £ [1,4] and 
equal-mass binaries with % = 5,/M 2 £ [—0.5,0.5]. We conclude that for the cases we study simulations must 
contain between three (in the equal-mass nonspinning case) and ten (the % = 0.5 case) orbits before merger, but 
there is also evidence that these are the regions of parameter space for which the least number of cycles will be 
needed. 

PACS numbers: 



I. INTRODUCTION 

Numerical simulations play a key role in efforts to de- 
tect gravitational waves (GWs) from compact-binary coales- 
cences, in particular from black-hole binaries. One of their 
most direct applications for GW searches has been as input 
in phenomenological [1-6| and effective-one-body (EOB) [7 - 
[T5l waveform models, which can in turn be used to con- 
struct template banks of theoretical waveforms for use in 
matched-filter searches in detector data. As part of the NINJA 
project [16 1 they are also used to test a battery of search 
pipelines: the numerical-relativity (NR) waveforms are in- 
jected into simulated detector noise and the search pipelines 
attempt to find them lfT71[T8l . 

Waveform models are constructed by combining informa- 
tion from post-Newtonian (PN) calculations of the long slow 
inspiral of the binary 1 19], and NR simulations of the last or- 
bits and merger. The PN approximation becomes less accu- 
rate as the binary approaches merger, and we will argue in 
this paper that at the current state of the art, PN errors domi- 
nate significantly over errors in the NR waveforms. Therefore 
more accurate waveform models can be produced by using 
less PN inspiral cycles, and correspondingly more NR cycles. 
The question we begin to address in this paper is: how many 
NR cycles are necessary in order to produce waveform mod- 
els that are sufficiently accurate for GW detection? By this 
we mean that the mismatch error in the full waveform (mini- 
mized with respect to all search parameters) is less than 3%, 
as described in Sec.llll 

This question is important because NR waveforms from 
an extremely large number of binary configurations may be 
necessary to produce waveform models that accurately repre- 
sent the entire black-hole-binary parameter space. NR simu- 
lations are computationally expensive, and the computational 



cost grows drastically with the length of the simulation. If the 
number of inspiral cycles in a simulation is to be increased by 
a factor of two, for example, the simulation doesn't only have 
to run for twice as long. To accurately capture the phase evo- 
lution of the binary over this extended time, higher numerical 
resolutions are required, and the overall increase in the cost of 
the simulation may be by factors of ten, both in memory usage 
and the time the simulation takes to run. 

Current "long" simulations cover ~10 GW cycles before 
merger |j20| , and the longest simulation to date covers about 
30 cycles |2T1 . The task of exploring the full black-hole- 
binary parameter space with long numerical simulations is 
a large-scale computational challenge that will take several 
years with current methods. On such a time scale, it is crucial 
to understand as much as possible about any requirements that 
may affect the overall cost by an order of magnitude! 

The necessary length of NR waveforms will change with 
improvements in the accuracy of PN methods and increased 
understanding of EOB methods. As such, we can only make 
a first step in considering waveform length requirements. Our 
method is to study hybrid waveforms, which are a simple con- 
nection of PN and NR results, and allow us to focus on the 
effects of the PN errors without having to take into account 
any additional artifacts that might arise in a given waveform- 
model construction procedure. We believe that this allows us 
to make the most conclusive statements that are possible with 
the current state of the art. 

We summarize mismatch calculations in Sec. [TTJ describe 
our numerical waveforms in Sec. Ill and the production of 
hybrid waveforms in Sec. |IV] We use hybrids based on two 
time-domain PN approximants, the TaylorTl and TaylorT4 
approximants, the accuracy of which we have already studied 
over a 10-cycle range before merger in [22 1. We put our pre- 
vious PN-NR comparisons into the context of GW searches 
by computing the mismatch between hybrids with differing 
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numbers of TaylorTl and NR cycles in Sec. |Vj and demon- 
strate that the hybrid errors are indeed dominated by PN ef- 
fects. We then turn to the question of length requirements for 
NR waveforms. 



We will explain and justify our approach in Sec. VI Since 
it has already been shown that pure PN waveforms truncated 
before merger are sufficient to detect nonspinning binaries 
with masses up to 12M [123 ], we consider binaries with to- 
tal mass M > 1OM , and our analysis is with respect to the 
Advanced LIGO detector [24]. Our study covers a significant 
portion of the parameter space that has already been treated 
by NR simulations: nonspinning binaries with mass ratios up 
to q = M2/M1 = 4, and equal -mass binaries with equal, non- 
precessing spins up to % = ±0.5. Results for nonspinning bi- 



naries are given in Sec. VII and spinning binaries are treated 



in Sec. VIII We will not consider phenomenological or EOB 
methods, and we will only briefly touch on parameter estima- 



tion in Sec. IX Our first conclusions about waveform length 
requirements are summarized and discussed in the Conclu- 
sion. 



II. MISMATCHES AND GW DETECTION 

The quantity that we will use to assess the detectability of 
our hybrid waveforms is the mismatch. The standard require- 
ment for current GW searches is that the members of a tem- 
plate bank be separated such that the no more than 10% of 
signals would be lost due to the mismatch between the signal 
and the template, and this translates into a mismatch between 
neighboring templates of less than 3% 11251 l26l . However, 
one must also take into account the mismatch due to errors 
in the theoretical waveforms; ultimately we want the sum of 
the template-spacing mismatch and the waveform-error mis- 
match to be less than 3%. We are then free to decide how to 
divide that mismatch between template spacing and waveform 
error. For example, the authors of the waveform accuracy 
study in [27] request that the waveform-error mismatch be be- 
low 0.5%. We could alternatively argue that the waveform 
error can be close to 3%, and the template spacing should be 
made correspondingly small. This question deserves further 
attention, but in this work we will consider 3% as the basic re- 
quirement, and also quote results consistent with a maximum 
mismatch requirement of 1.5% and 0.5%. 

For two GW signals hi and hi, we define an inner product 
in the Fourier domain weighted by the power spectral density 
of the detector noise, S„(f), as 
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where [/ m i n ,/ m ax] is the intersection of the chosen sensitivity 
range of the detector ([20, 10 4 ] Hz in this study) and the range 
of validity of the waveform data (in most cases fM € [1 .25 x 
10~ 3 ,0.15]). hi(f) and hi{f) denote the Fourier transform 
of h\ and h 2 , respectively. Our data in general represent the 
Weyl scalar ^(f), not the wave strain h(t), but the two are 
related by ¥4 = h + — ih x , and two time integrations can be 
performed trivially in the frequency domain. 



Given the definition of the inner product (h\\h.2), we nor- 
malize it and maximize over phase and time offsets in the 
data. This is the faithfulness of the waveform: it is a mea- 
sure of how "far" a theoretical waveform is from a suppos- 
edly true waveform with the same physical parameters. The 
faithfulness mismatch is the deviation of the faithfulness from 
unity: 



1 — max — , 



(2) 



In a GW search, the goal is essentially to find the template- 
bank member that has the smallest mismatch with the detec- 
tor data. We are therefore really interested in the mismatch 
optimized over all of the members of a theoretical waveform 
family, i.e., minimized with respect to the binary's intrinsic 
physical parameters (the mass, mass ratio and spins), and the 
extrinsic parameters (the position of the binary in the sky and 
its orientation). The fully optimized mismatch is the same as 
(1 —F), where F is the waveform's fitting factor. In our study 
we have access to signals of isolated binary configurations, 
and the only parameter we can optimize with respect to is the 
total mass. As such, an optimization of the mismatch with re- 
spect to total mass can be no more than an upper bound on the 
full optimized mismatch. But we will argue in Sec. | VII| that 
optimization with respect to the other physical parameters will 
not qualitatively alter our results. In Sec. IV and[V] where we 
assess the sources of physical error in our hybrids, we will 
not optimize with respect to total mass. 



VI 



But in Sees. 

|IX| where we make calculations relevant to GW detection, we 
will use the mass optimization, and intend these results to be a 
reasonable estimate of the mismatch calculated from the fully 
optimized fitting factor. 

Our main focus in this work is GW detection, but we will 
also refer to a quantity that is relevant to parameter estima- 
tion. If we have two theoretical waveforms for the same phys- 
ical system, h\ and h 2 , each with some associated uncertainty, 
then we can define Sh = h\ — h 2 and then ||5/i|| 2 = (Sh\Sh). 
If 1 1 8h 1 1 < 1, then the two waveforms would be indistinguish- 
able in a detector measurement |27|. The indistinguishability 
of two waveforms depends on the signal-to-noise ratio (SNR) 
of the detection: if the SNR is sufficiently low, then any two 
waveforms will be indistinguishable (although below an SNR 
of eight they will not be considered detectable anyway), and 
for two waveforms of arbitrarily low we can always 

find an SNR high enough such that they can in fact be distin- 
guished. For example, in 1 29 1 we found that current NR equal- 
mass nonspinning waveforms from five different codes would 
be indistinguishable in the Advanced ground-based detectors 
if the SNR is below 25 . And in [ 5 ] it was shown that q = 2 hy- 
brid waveforms produced with the BAM and Llama ll30l [3~TTl 
codes are indistinguishable for an SNR below roughly 20. 

The mismatch error between two waveforms can be easily 
related to \\8h\\ |32|. If p = \\h\\ is the optimal SNR, then 
1 1 8h 1 1 /p 2 « 2JK. This means that if two waveforms meet 
the detection criteria of ^0 < 0.03, then they will be indistin- 
guishable for p < 4, which is too weak a signal to be detected. 
If we want the waveforms to be indistinguishable at an SNR 
of p > 8, then the mismatch must be below about 0.8%. We 
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see, then, that the accuracy requirements for two waveforms 
to be indistinguishable are in general far more stringent than 



those for detection. We will return to this point in Sec. IX 



All of the results in this work will be with respect to the 
Advanced LIGO detector [24 33], and in general we will use 
a low -frequency cut-off of 20 Hz. 



III. NUMERICAL WAVEFORMS 

We consider two families of black-hole binaries. The first is 
equal-mass binaries in which the spin of each black hole is the 
same, Si jM\ = S^jM^ — %, and the spins are parallel or anti- 
parallel to the orbital angular momentum. In these configura- 
tions the spins do not precess, making this a simple subfam- 
ily of the black-hole-binary parameter space. The spins used 
were % = {0, ±0.25, ±0.50, ±0.75, ±0.85}, although we will 
only provide length requirements up to % = ±0.5 for reasons 



that will be explained in Sec. VIII The second family is non- 
spinning binaries with mass ratios q =M2/M\ = {1,2,3,4}, 
where we labelled the individual masses such that Mi > M\ . 

The simulations were produced using the BAM code IT341 
[35) . They cover 6-10 orbits before merger. The phase error 
of the waveforms is at most 0.15 rad during the inspiral phase 
(up to Mco = 0.1), and on the order of 1 rad during merger and 
ringdown. The amplitude accuracy is within 1% during the 
inspiral, and within 5% during the merger and ringdown. Of 
more relevance to GW detection is the mismatch error of the 
waveforms, which is below 10~ 4 for all cases. Note that this 
is well within the detection requirement of 0.03, and indeed 
these waveforms are well within the accuracy requirements 
for both detection and parameter estimation with current and 
planned ground-based detectors ll29l . These waveforms were 
all presented in detail in | 



IV. HYBRID INSPIRAL-MERGER-RINGDOWN 
WAVEFORMS 

Complete inspiral-merger-ringdown waveforms that in- 
clude arbitrarily large numbers of GW cycles before merger 
can be constructed by connecting PN and NR results fl~l|2][5] 
[38l l39l . We model the inspiral regime using the time-domain 
approximants TaylorTl and TaylorT4 |5][37). For nonspin- 
ning binaries the phase evolution is described to 3.5PN order, 
and the amplitude is given to 3PN order fiOl . For spinning 
binaries, spin effects in the phase evolution are included only 
up to 2.5PN order (37), and in the case of TaylorT4 we adopt 
two approaches to truncating the expansion in a consistent 
way f5] [22); the amplitude includes spin contributions up to 
2PN order ED. 

The PN results are expansions in the frequency, x = 
(MH) 2 / 3 , where £2 is the orbital frequency of the binary mo- 
tion. When decomposing the GW signal in spherical harmon- 
ics, £2 is related to the frequency co of the (I = 2,m = ±2) 
modes by |o| = 2D.. As such, the PN waveforms are most 
accurate for small x, i.e., many orbits before merger. As the 
binary approaches merger, the errors in the PN approximation 



grow. The convergence properties of the PN expansion are not 
fully understood, and it is not possible to provide clear error 
estimates. One way to gain some insight into the PN errors 
is to compare results from different PN approximants and at 
different PN orders, and we will return to this idea later. 

The only way to definitively quantify the PN errors is to 
compare with fully general relativistic results. This can be 
done over the (relatively small) number of GW cycles for 
which we have both PN and NR waveforms. PN-NR com- 
parisons of phase and amplitude have been performed over 
~ 10 orbits prior to merger OTl [36l I42H441 . Comparisons for 
the cases we consider here are given in [22]. There we found 
that the TaylorT4 approximant has an accumulated phase er- 
ror of no more than 0.2 rad for nonspinning cases, but the 
phase error grows to up to 2 rad for equal-mass cases with 
large spin, in particular the % < cases. The behaviour of 
the TaylorTl phase appears to be most consistent across the 
parameter space that we consider, with an accumulated phase 
error of about 1 rad in all cases, although in most cases a ver- 
sion of TaylorT4 is more accurate. The PN amplitude (at the 
highest-known PN order) has an error as high as 4% in the 
X = 0.85 case, and drops to around 2% in the % = —0.85 
case, while for nonspinning configurations it is around 3%. 
(Note however that the uncertainty in the NR amplitude for 
all of these cases is 1%.) The existence of this region of good 
agreement between PN and NR results is what allows us to 
combine the two into hybrid waveforms. 

We construct time-domain hybrid waveforms by first de- 
composing the PN and NR waveforms into their amplitude 
[Apn(?) andANR(f)] and phase [0pn(O and 0nr(?)]. From the 
phase we can in turn calculate the frequency, co(t) = d t (j)(t). 
We then choose a matching frequency, co m , and determine the 
time when both 0pn and 0nr reach that frequency, and con- 
nect the two phase functions at that time. We independently 
determine the time when the PN and NR amplitudes agree, 
and connect the amplitudes at that time. In this procedure 
we do not require that the transition between the PN and NR 
phases and amplitudes is any more than continuous. This pro- 
cedure is performed on ^4. In mismatch calculations we first 
calculate ^4 by making a fast Fourier transform (FFT) of ¥4, 
and then produce h by dividing by -co 2 . 



A. Errors due to the hybridization process 

There are currently several procedures to produce hybrid 
waveforms that have been used in the literature, based either 
on a matching in the time or Fourier domains |fTll2"ll5l l38l[39l . 
Each procedure will itself introduce artifacts into the wave- 
form, and we may be concerned that these artifacts will be 
a large source of error. Here we will compare hybrids pro- 
duced using three different methods, and show that in fact the 
mismatch error introduced by the hybridization procedure is 
negligible. 

We consider three equal-mass nonspinning hybrids. The 
first is produced by the method that we have just described. 
The second is produced using the method described in JT] [2], 
but applied to ^4 instead of the wave strain h, because we 
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FIG. 1 : Mismatch between equal-mass nonspinning T4+NR hybrids 
constructed using two different time-domain methods (solid line), 
and between hybrids constructed in the time domain (TD) and the 
frequency domain (FD). The time-domain hybrids are produced with 
matching frequency MOm = 0.07, and the frequency-domain hybrid 
is matched at Mco = 0.079. 



want to assess the errors due to the hybridization process, 
without any contamination by errors that may be introduced 
by a time-domain integration of VP4 to h. Our specific method 
is to choose a matching frequency, ©„,, and to then locate the 
time in both the PN and NR waveforms when that frequency is 
reached. We then combine the two waveforms over a 200M- 
long window, aligning the waveforms such that the quantity 

A¥ = jf* (v 4 Mt)-ae i5 ^4Mt + $t)y ' dt (3) 

is minimized, where ?o and t\ are respectively 100M before 
and after the time t m at which each waveform reaches co m , a is 
a scale factor, 8t and 50 are time and phase offsets, and the 
waveforms are initially aligned so that t m is the same for both. 
The hybrid is constructed by making a linear transition be- 
tween ^.pn and ^.nr over the matching window. For both 
hybrids, we choose a matching frequency of Ma>,„ = 0.07. 

The third hybrid is constructed in the frequency domain, 
using a variant of the method described in 0. We produce 
an FFT of the time-domain TaylorT4 approximant, and an 
FFT of the numerical ¥4 data. The phase of the frequency- 
domain PN and NR signals is then matched in the window 
Mco £ [0.0566,0.113] and the continuous transition is carried 
out at the matching frequency Mco m = 0.079. 

The non-mass-optimized mismatch is shown in Fig. [T] We 
see that the maximum mismatch is about 0.025% between 
the two time-domain hybrids, and 0.03% between the time- 
domain and frequency-domain hybrids. Clearly the error due 
purely to the hybridization procedure is, like the mismatch er- 
ror of the numerical waveforms, negligible. 

If we consider the indistinguishability criteria, \\8h\\ < 1, 
then these hybrids would be indistinguishable for SNRs of 
p < 40. 

We note, however, that the difference between hybrids 
constructed with ostensibly the same numerical waveforms 
and PN approximants may have larger differences than those 



shown here. For example, if we compare either of the hybrids 
we have just described, with hybrids constructed using the in- 
tegrated wave strain, as in 1 1, 2 1, then the mismatch can be as 
high as 0.8%. This is due not to the hybridization process, but 
to artifacts introduced in the time-domain integration of ¥4 to 
h. 



V. COMPARISON BETWEEN PN AND NR DURING THE 
LATE INSPIRAL 

Direct comparisons between the phase and amplitude of PN 
approximants and NR waveforms have been made in I12T1 [36l 
I42ti46l . and we refer the reader to ll22l for PN-NR compar- 
isons of the binary configurations studied here. In this section 
we reframe those comparisons in the context of hybrid wave- 
forms and mismatches. 

For several choices of binary configuration, we construct 
a fiducial reference hybrid waveform by matching a 3.5PN 
TaylorTl waveform to the longest available NR simulation. 
We then construct a set of "candidate" hybrid waveforms by 
sliding the matching region over the NR waveform (such that 
progressively shorter NR segments are employed in the hybrid 
construction). How well one candidate hybrid waveform can 
be distinguished from the reference waveform is given by the 
mismatch between the "reference" and the "candidate". 

We study three cases, employing three different NR sim- 
ulations: an equal-mass, non-spinning simulation performed 
using the SpEC code (47] 08), and two of the BAM simu- 
lations presented in 11221 . equal-mass simulations with spins 
aligned and anti-aligned to the orbital angular momentum 
(X = ±0.75). Note that in this section, we choose a low- 
frequency cut-off of 10 Hz, the mismatches are not maximized 
over the binary mass, and the hybrids are constructed using the 
procedure described in flT] [2] . 

Figure [2] shows the mismatch between the reference hybrid 
waveform and the hybrid waveforms constructed using differ- 
ent matching regions in the equal-mass, non-spinning case. 

The matching frequency of the reference waveform is 
McOyn.r = 0.031 - 0.038. We will denote the matching fre- 
quency of the candidate hybrid by (Om. 

When <a m r = (Mm, then the two waveforms are identical, and 
the mismatch is zero. As the matching frequency of the can- 
didate waveform 0) m is increased, the mismatches grow. This 
gives us a picture of the mismatch due to the disagreement of 
TaylorTl and full numerical results over the frequency range 
Mco £ [tOm.r, 0) m ]. For the maximum matching frequency we 
choose, Ma^ = 0.075, the maximum mismatch is 3.5%. 

Figures[5]and[4]show the mismatches between the reference 
and candidate hybrid waveforms for the case of equal-mass bi- 
naries with aligned and anti-aligned spins % = ±0.75, respec- 
tively. The x = —0.75 reference waveform was constructed 
with a matching frequency of Mo m r ~ 0.038 — 0.044, and the 
X = 0.75 waveform is matched at Mco mj ~ 0.057 - 0.063. 

We see from these results that the mismatch rises above 
3% at lower matching frequencies than in the non-spinning 
example. This is a reflection of the fact that the TaylorTl 
approximant performs worse in the spinning cases. 
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FIG. 2: Mismatch between the reference hybrid waveform and the 
hybrid waveforms constructed using different matching frequencies. 
The hybrid waveforms are constructed by matching TaylorTl PN 
waveforms with the Caltech-Cornell equal-mass, non-spinning NR 
simulation. 



FIG. 4: Same as in Fig.[2]except that the hybrid waveforms are con- 
structed by matching TaylorTl waveforms with BAM simulation of 
equal-mass binary with spins % — —0.75. 



VI. WAVEFORM LENGTH REQUIREMENTS: GENERAL 
PROCEDURE 
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FIG. 3: Same as in Fig.|2]except that the hybrid waveforms are con- 
structed by matching TaylorTl waveforms with BAM simulation of 
equal-mass binary with spins % = 0.75. 



These first results give us an indication of the mismatch er- 
ror associated with the disagreement between the TaylorTl 
approximant and NR results over the frequency range for 
which NR results exist. They show that the accumulation of 
mismatch error over that frequency range due to PN errors 
is orders of magnitude larger than the mismatch error due to 
errors in the numerical waveforms or to the hybridization pro- 
cess. These mismatches could drop significantly if we were 
to minimize the mismatch with respect to binary mass, but 
that would not alter the qualitative observation that PN errors 
dominate the error in our hybrid waveforms. 



In producing hybrid waveforms, we have seen that the dom- 
inant error source is the error in the PN approximation. We 
expect the PN waveforms to be increasingly accurate as we 
move to lower frequencies (i.e., earlier in the inspiral), and 
so the lower the frequency co m at which we can match (i.e., 
the more NR cycles we can use) the better. Our goal now is 
to determine the highest acceptable matching frequency for a 
variety of accuracy criteria for the full waveforms. 

A previous study showed that pure PN approximants of 
nonspinning binaries (i.e., without any stitching to merger 
and ringdown signals), are adequate for detection purposes 
up to about 12M (23l. Therefore our PN+NR hybrids have 
to be accurate only for total masses M > 1OM . Above 
these masses the mismatch between PN inspiral waveforms 
and full inspiral-merger-ringdown waveforms deteriorates 
quickly EH M ED . 

The main difficulty in assessing the physical accuracy of 
hybrids is that, since we cannot rigorously quantify error bars 
in the PN waveforms, any error estimates we make on our 
complete waveforms can be no better than an educated guess. 
The best we can do is to make clear what assumptions have 
gone into our guess, and to what extent this guess can be con- 
sidered conservative or optimistic. We will now describe the 
procedure we have used, and our justifications for it. 

Our approach is the following. We start by defining a model 
waveform that we will regard for the purposes of this study as 
the "true" physical waveform. 

For nonspinning binaries we know that during the fre- 
quency range in which PN and NR results overlap, the Tay- 
lorT4 approximant captures the phase evolution with high ac- 
curacy. We now assume that the TaylorT4 phase continues 
to be physically correct to much lower frequencies, and we 
construct a TaylorT4+NR hybrid and treat that as the true 
GW waveform from our system. In other words, we treat 
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this hybrid as if it were a numerical waveform starting at 
Mco = 0.006, which corresponds to 20 Hz for a 10M Q binary. 
Such a waveform covers ^650 cycles before merger. We have 
no expectation of ever producing such a waveform in a numer- 
ical simulation, but for the purposes of this exercise we will 
treat our hybrid as if it were the result of just such a simula- 
tion, and refer to it as NRl- 

Our goal is to assess the accuracy of PN+NR hybrids con- 
structed using a reasonable choice of PN approximant. Based 
on its phase accuracy near merger, we choose the TaylorTl 
approximant. We construct hybrids of TaylorTl with NRl, 
and due to the great length of NRl we have the freedom of 
matching at arbitrarily low frequencies. 

We then calculate the mismatch between T1+NRl and the 
true NRl waveform, and use this to estimate the ability of 
T1+NRl to detect the true signal. This process is repeated for 
a range of matching frequencies, which allows us to determine 
the highest matching frequency at which the T1+NRl hybrid 
is sufficiently accurate for detection of binaries above 10 M e . 

The one key assumption in this procedure is that the Tay- 
lorTl phase can be trusted to much lower frequencies than 
those where it has currently been compared with full numer- 
ical results. The first and most detailed study of TaylorT4 
was presented in iPTI . which included frequencies down to 
Mco ps 0.035, so we know that TaylorT4 (at least in the non- 
spinning case) is accurate down to that frequency. But below 
that frequency we have no information about its accuracy with 
respect to waveforms from full general relativity. We expect 
that the difference between the TaylorT4 and TaylorTl phases 
provides a reasonable estimate of the phase error in a typical 
approximant, but there is no way to prove this. Ultimately 
one is reduced to a statement of faith in the accuracy of PN 
methods. 

Let us illustrate this point with two extreme views. 

Instead of comparing TaylorTl and TaylorT4 hybrids, we 
could compare hybrids produced using different effective- 
one-body (EOB) calculations. It has been claimed that EOB 
methods, particularly after calibration to numerical results, 
provide an extremely accurate model of the full GW signal J7J- 
[131 15014541 . As with standard PN appro ximants, EOB results 
can only be tested in the regime where NR results exist. How- 
ever, all PN and EOB results make use of expressions for the 
GW flux and energy loss of the binary, and it has been shown 
in 1 55 1 that EOB estimates of these quantities agree with NR 
results far better than any standard PN estimates. We might 
then expect, on this evidence alone, that EOB results are far 
more accurate than any standard PN approximant, and the ap- 
propriate comparison would be between hybrids produced us- 
ing variants of the EOB method. 

This view suggests that the comparison we are proposing 
— between Tl and T4 hybrids — is overly pessimistic of the 
potential physical fidelity of hybrid waveforms. 

An alternative view is that EOB waveforms (whether the 
final waveform phase, or the flux and energy-derivative ingre- 
dients) have nonetheless only been compared with full GR re- 
sults close to merger, and their accuracy at lower frequencies 
is unknown. One could argue that we should be much more 
conservative, and instead of comparing 3.5PN TaylorTl and 
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FIG. 5: Mismatch between T4+NR hybrids constructed with match- 
ing frequencies Mco m = 0.06 and M(O m = 0.09. We see that the mis- 
match between T4 hybrids constructed at different matching frequen- 
cies is below 0.12% over the entire mass range that we consider. 



TaylorT4 hybrids, we should compare hybrids constructed us- 
ing different PN orders. Only by comparing, for example, 
2.5PN and 3.5PN results, can we hope to estimate the error 
between 3.5PN results and those from full GR. As we will see 



in Sec. VII this provides an extremely pessimistic view of the 
accuracy of 3.5PN waveforms. 

In the end, then, we consider our proposed T1-T4 compari- 
son to be a reasonable compromise between these two extreme 
views. But we hope this long preamble serves as a caveat that 
all such comparisons are by no means fully conclusive. 



VII. NONSPINNING BINARIES 

We will now investigate the length requirements for simu- 
lations of nonspinning binaries. We start with the equal-mass 
nonspinning case. 

We first construct a hybrid between 3.5PN TaylorT4 and 
our numerical waveform. We make the transition between PN 
and NR at a matching frequency Mco m = 0.09. Fig. [5] shows 
the mismatch between this waveform and a hybrid produced 
with a different matching frequency, and suggests that for our 
purposes the resulting hybrid does not depend strongly on the 
matching frequency (O m , since we are not concerned by mis- 
matches that are below 0.5%. Note that these results are a 
further reflection of the good agreement of TaylorT4 with the 
numerical phase at these frequencies; we see from Fig.|2]that 
the same would not be true for Tl+NR hybrids. 

We treat this waveform as if it were an extremely long nu- 
merical waveform, and denote it NRl- This is our target GW 
signal. It is around 4 x 10 5 M in duration, making it two or- 
ders of magnitude longer than the longest current equal-mass 
nonspinning waveform ETl . TaylorTl+NRL hybrids are then 
constructed with a range of matching frequencies, and com- 
pared with NRl, as described in the previous section. 

Fig. [6] shows the mismatch between NRl and T1+NRl 
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FIG. 6: Mismatches between T1+NRl and NRl hybrids, for match- 
ing frequencies Mco m = {0.045,0.05,0.06,0.07,0.09}. 



FIG. 7: Mismatches between T1+NRl and NRl hybrids, where the 
TaylorTl approximant is evaluated at only 2.5PN order. 



hybrids for Mco m = 0.045,0.05,0.06,0.07,0.09. Above the 
highest matching frequency of Mco m = 0.09 both hybrids 
are identical for all the matching frequencies we have cho- 
sen. This means that for masses above about 150M Q , the 
two waveforms will be almost identical within the sensitiv- 
ity range of the detector, and their mismatch should approach 
zero. For this reason Fig. [6] considers masses only up to 
1OOM , and indeed we see that the mismatch has essentially 
dropped to zero by 100M Q . 

At lower masses, the difference between the two PN ap- 
proximants dominates, and the mismatch rapidly grows. At 
yet lower masses, the PN approximants are being sampled at 
lower frequencies, where their agreement is better, and we ex- 
pect the mismatch to fall. This behavior can be seen for the 
matching frequency Mco m = 0.09, where the peak mismatch 
is at about 13 M . For lower matching frequencies, the peak 
mismatch is pushed to lower masses, and for Mco m < 0.05 we 
do not see it in the figure. We expect that the location of the 
mismatch maximum is related to the mass at which Mco m is at 
the detector's most sensitive frequencies. The detector's peak 
sensitivity is at around 200Hz, and for example Mco = 0.09 
corresponds to that frequency for M = 15.4 M , which is con- 
sistent with Fig. [6] 

As described in Sec. |VI| the mismatch is optimized with 
respect to the total mass of the binary. We cannot optimize 
over the mass ratio, because we do not have access to numer- 
ical waveforms with q close to unity. However, we can make 
a crude estimate of the magnitude of the effect of mass-ratio 
optimization. If we consider a 1OM binary, then most of the 
hybrid in the detector's sensitivity band is from the PN contri- 
bution. We vary the mass ratio of a TaylorTl waveform until 
its phase agreement with TaylorT4 is optimized, and then con- 
struct a TaylorTl+NPvL hybrid using that mass ratio, and re- 
peat the above analysis with Mco m = 0.06. The optimal mass 
ratio is q = 0.99985, and the change in the mismatch is indis- 
tinguishable on the scale of Fig. [6] This suggests that the mass 
optimization by itself is sufficient to give an indication of the 
full mismatch error of our hybrids. 

If we are willing to consider waveforms acceptable so long 



as the waveform-error mismatch is below 3%, then we see 
from Fig.|6]that the Mco,„ = 0.09 Tl hybrids are almost good 
enough, and hybrids produced with any lower matching fre- 
quency are acceptable. With this criterion, we need wave- 
forms that contain only five GW cycles before merger, or 
about two orbits! (Note that the number of GW cycles before 
the peak amplitude is more than simply twice the number of 
orbits, because the approximate quadrupole relation co = 2Q, 
no longer holds.) 

Being more realistic, and demanding that the mismatch is 
below 1.5% (so that the remainder of our allowed 3% mis- 
match is taken up by the template spacing), a matching fre- 
quency of Mco m = 0.06 is more than acceptable. This means 
that we need about 11 GW cycles before merger, or about 
five orbits. If we are especially stringent and require that the 
error mismatch be below 0.5%, then we must match below 
Mco,„ = 0.05, and this requires 15 cycles (seven orbits). 

For the equal-mass nonspinning case, then, we see that the 
NR waveform length requirements for GW detection are not 
very high: five orbits are sufficient, and seven orbits are plenty. 

At this point we may ask what happens if we decide to be 
much more conservative in estimating the PN error. We can 
do this by repeating the exercise above, but instead using Tay- 
lorTl waveforms at only 2.5PN order. The results are shown 
in Fig. [7] We see that the variation in the mismatch is dra- 
matic. Even if we choose a matching frequency of Mco,„ = 
0.04 (which corresponds to 23 GW cycles before merger), the 
waveforms are not accurate enough for searches, irrespective 
of template bank spacing, for masses below 15 M©. Only if we 
match at Mco = 0.03 (40 cycles before merger) are the wave- 
forms usable, and to get close to the more rigorous mismatch 
requirement of 1.5% for all masses, we need Mco = 0.02, 
which corresponds to 80 GW cycles before merger. If one is 
to adopt the view that this is a realistic estimate of the uncer- 
tainty in our 3.5PN approximants, then we must conclude that 
hybrids constructed from even the longest numerical wave- 
forms currently published are only useful for searches down 
to 15-20 solar masses. It is the opinion of the authors, how- 
ever, that this is a gross exaggeration of the error in 3.5PN 
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TABLE I: Empirically calculated 3PN and 3.5PN coefficients that 
produce a TaylorT4 approximant that agrees well with the NR phase 
over the last ten cycles up to Mco — 0. 1 . 



X 


3PN 


3.5PN 


-0.85 


989.4 


16.26 


-0.75 


775.3 


21.11 


-0.50 


259.2 


23.99 


-0.25 


39.73 


-1.22 





0.868 


-0.111 


0.25 


-128.9 


0.320 


0.50 


-204.6 


1.070 


0.75 


-298.5 


38.97 


0.85 


-268.6 


-1.009 



VIII. EQUAL-MASS SPINNING BINARIES 



approximants ! 

We now consider higher mass ratios. The results in 
suggest that TaylorT4 is also a good model for the true wave- 
form for the mass ratios that we consider, q = {1,2,3,4}. We 
find that the maximum acceptable matching frequency drops 
as the mass ratio increases, so that we need longer numerical 
waveforms. Fig. [8] shows the results for q = 4, which is the 
most extreme case. The figure shows results with matching 
frequencies Mco m = {0.04,0.045,0.05,0.06,0.07}. Clearly 
matching frequencies of Mco — 0.07 will not be sufficient even 
if we allow mismatches up to 3%, butMo,,, = 0.06 (15 cycles, 
or about 7 orbits before merger) is borderline. A matching 
frequency of Mco,„ = 0.05 (21 cycles, 10 orbits) ensures mis- 
matches no higher than 1.5%, and mismatches below 0.5% 
require a matching frequency of Mco — 0.04 (33 cycles, or 
about 15 orbits). 

We see, then, that for higher mass ratios it is not suffi- 
cient to produce waveforms of the same modest lengths as 
in the equal-mass case. This is unfortunate, because higher- 
mass-ratio simulations are far more computationally expen- 
sive. This result also highlights a drawback of direct com- 
parisons between NR and PN results over a small number of 
cycles before merger: in such comparisons, the performance 
of TaylorTl and TaylorT4 is roughly the same for q = 1 and 
q = 4, and we might therefore conclude that the length re- 
quirements for NR waveforms would be the same if we want 
to produce sufficiently accurate TaylorTl+NR hybrids in ei- 
ther case. But these results show otherwise; we need more 
cycles for higher mass ratios. 

The q = 4 numerical waveform used for this study covers 17 
GW cycles before merger, and so would in principle be usable 
for searches down to 10M R , although a slightly longer wave- 
form (for example including 10 clean inspiral cycles) would 
be preferable. It is clear, however, that (1) producing accept- 
able q = 4 waveforms is certainly feasible with current codes, 
but (2) it is non-trivial to estimate the length requirements for 
yet higher mass ratios. 



We now tackle the cases with spin. The TaylorT4 approx- 
imant no longer accurately tracks the true phase in the last 
cycles before merger, and it would be difficult to justify its 
use in constructing an NRl hybrid as in the previous section. 
(See Figs. 7 and 8 in ED.) 

Without a promising approximant to use to produce a stand- 
in for a true waveform, our next-best option is to construct 
one. The nonspinning terms of TaylorT4 are known up to 
3.5PN order, but the spinning terms are known only up to 
2.5PN. Perhaps if we knew the 3PN and 3.5PN contributions, 
the spinning TaylorT4 would be just as good as the nonspin- 
ning version? Whether this is true or not, we can certainly 
introduce 3PN and 3.5PN terms that improve its agreement 
with numerical results. Our procedure is to adjust the coef- 
ficients of such terms so as to minimize the square-integral 
phase difference between the PN and NR waveforms, as de- 
fined by | 



A<j)(t N ) = 



1 



, s2 • 

0nr(O-0pn(O dt 



',v 



1/2 



(4) 



where N is the number of cycles included in the comparison, 
and we choose N = 10. We find in practice that the phase 
evolution is almost identical using an approximant with either 
an optimized 3PN coefficient only, or both 3PN and 3.5PN 
coefficients. This suggests that our fitting procedure over 10 
GW cycles is not very sensitive to higher order coefficients. 

However, the minimization process is sensitive to the initial 
guess for the coefficients, and it is possible that other choices 
are possible with similar results. The particular coefficients 
that we use are given in Tab. [I] We note that the 3PN coeffi- 
cients depend roughly monotonically on the spin of the black 
holes, but we will not attempt to infer any significance on 
the particular values that we obtain. It is important to bear 
in mind also that these coefficients are chosen to achieve the 
best agreement over only ten GW cycles, and those cycles are 
near merger, where the PN approximation is close to break- 
ing down. The "correct" 3PN and 3.5PN coefficients, when 
they are derived analytically, will be such that they will be 
expected to lead to a GW phase that is physically correct at 
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low, not high, frequencies. The coefficients that we are using 
may have the reverse properties, and lead to a poor estimate 
of the phase at low frequencies — and in fact our ad-hoc co- 
efficients may adversely distort the phase function during the 
early inspiral. If anything, this will lead us to conclude that 
our waveforms should be much longer than they really need 
to be, and so our results will still provide an upper bound on 
the necessary waveform length. And this is our goal. 

This point is illustrated in Fig. |9| which shows the accumu- 
lated phase disagreement between the optimized 3.5PN Tay- 
lorT4 approximant, and our three standard approximants, Tay- 
lorTl, TaylorT4, and TaylorT4 tluncate d for the % = 0.5 case. 
The results in ll22l show that T4 agrees best with the NR 
data over the ten cycles of comparison, Tl is the next best, 
and T4 trun cated performs worst. But if we look now at Fig. [9] 
which shows the phase disagreement during all of the early 
inspiral up to Mco = 0.06, we get a different impression. Now 
the T4 truncate( j approximant performs best, while T4 is next 
best, and Tl is the worst. This confirms our suspicion that 
the relative performance of different approximants during the 
last orbits of the binary may be quite different to their per- 
formance during the earlier inspiral — but it also tells us that 
if we compare the hybrids produced by the optimized 3.5PN 
T4 approximant with TaylorTl hybrids, then this will give us 
the most conservative estimate of the length requirements of 
our waveforms. This is because it is the phase disagreement 
that dominates the mismatch, and Fig.[9]suggests that the Tl- 
T43 5pn mismatches will be the worst. 

We will focus on two cases, % = —0.5 and % = 0.5. The 



mismatch plots for these cases are shown in Fig. 10 The most 
notable aspect of these two plots is the dramatic difference in 
the mismatches as a function of matching frequency between 
the two cases. For the anti-hangup case % = —0.5, a matching 
frequency of about Mco m = 0.055 appears to be sufficient to 
achieve mismatches below 1.5%. This corresponds to about 
10 cycles before merger, and is comparable to what we saw 
in the nonspinning case. For the hangup case % = 0.5, on the 
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other hand, a matching frequency of between 0.04 and 0.05 is 
necessary. A matching frequency of Mco m = 0.04 corresponds 
to 28 cycles before merger, or 13 orbits. Some of this differ- 
ence is due to the effect of spin on the rate of inspiral: from 
a given frequency, there will be more cycles before merger 
in the hangup case than in the anti-hangup case. Since the 
hangup case requires matching at a lower frequency, this adds 
yet more cycles to our estimate — and so we end up with al- 
most a factor of three difference in the number of NR cycles 
that are required. 

To test the robustness of our results, we have also calculated 
the best-fit modifications to the 3PN and 3.5PN by optimizing 
over a smaller number of NR cycles. There is a large varia- 
tion in the parameters as the number of included NR cycles 
is changed, but we have found that the change in the overall 
mismatch results is at a level that would not be distinguishable 
in the figures. 

Tab. [II] summarizes the number of cycles required before 
merger for the different cases, for different levels of mismatch 
requirement. Spins up to only \%\ = 0.5 are used, because in 
higher-spin cases the deviation of the known T4 approximant 
from the NR phase is so large that the required 3PN and 3.5PN 
coefficients appear to distort the overall phase evolution too 
severely. We can infer from this table, however, that for % < 
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TABLE II: Summary of mismatch calculations. The second, third 
and fourth columns give the minimum number of numerical GW cy- 
cles before merger that ensure mismatches below 3%, 1.5% and 0.5% 
for all masses above 1OM . The last column indicates the lowest 
mass for which the numerical waveforms studied here could be used 
for searches, assuming that the mismatch error can be as high as 3%. 



Configuration 


Jt <3% 


Jt < 1.5% 


Jt < 0.5% 


M min /M Q 
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q = 3 
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q = 4 


15.0 
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—0.5, simulations that include 20 GW cycles before merger 
(10 orbits) should allow the construction of Tl hybrids that are 
within the 3% mismatch accuracy requirement, but for high 
X > 0.5 many more cycles may be required. 

The table also indicates, assuming that we will allow a max- 
imum error mismatch up to 3%, the minimum mass that can 
be searched for using hybrids constructed from the numerical 
waveforms presented here. These numbers are provided as a 
snapshot of what can be done with the waveforms that exist as 
of the writing of this paper. In most cases we can produce hy- 
brids that are acceptable for searches down to 1OM , but for 
some spinning cases we can use our hybrids to search down to 
only 15M Q . We hope that the results in this paper ultimately 
indicate upper bounds on the the required length of numerical 
waveforms for GW detections. It should be clear from Tab.HIl 
that an improvement in PN approximants (due for example to 
the calculation of higher-order spin terms) would have a sig- 
nificant effect on the NR waveform length requirements. If the 
accuracy of spinning approximants were improved to the level 
of their equal-mass nonspinning counterparts, then numerical 
waveforms covering only 5-7 orbits before merger would be 
sufficient for GW searches. 

On the other hand, we also see that the length require- 
ments vary significantly between binary configurations, and 
that more NR cycles are required as we approach the extremes 
of our parameter choices. We cannot determine from these 
results what the length requirements would be for q > 4 non- 
spinning waveforms, or even for highly spinning waveforms 
with q > 1 . In this sense, the "worst" cases remain to be stud- 
ied. 



As discussed in Sec. [TTJ if for our Tl and T4 hybrids we 
have | | < 1, then the two waveforms are indistinguishable. 
This means that the accuracy of the estimation of the intrin- 
sic parameters of the binary is determined by the SNR of the 
signal, and not by any error in the waveforms. In other words, 
if the waveforms are indistinguishable at the SNR of a given 
measurement, then the maximum parameter information can 
be extracted from that measurement, and is not limited by the 
accuracy of the waveforms. 

It is reasonable to expect SNRs as high as 30 in Advanced 
detectors, and for waveforms to be indistinguishable at that 
SNR, the mismatch error must be below 0.05% (see again the 
discussion in Sec. [II]). It is clear even in the best case (equal- 
mass nonspinning) that our hybrids meet this criteria only for 
binary masses higher than 20 M . To achieve such a low mis- 
match down to 1OM would require numerical waveforms far 
longer than any that have yet been produced. We have also 
seen that the mismatch in the hybrids due to artifacts from the 
hybridization process alone are at 0.03%, and so producing 
hybrids that are indistinguishable for parameter estimation up 
to an SNR of 30 is a challenge irrespective of the problems of 
PN errors and waveform length. 

However, even though the hybrids we have presented may 
not be indistinguishable at the potential SNRs of Advanced 
detectors, they can still be used to estimate the parameters of 
a binary, and the question remains what the accuracy of the 
parameter estimation can in principle be, and whether this ac- 
curacy is sufficient for GW astronomy applications in the near 
future. 

To give an indication of what these errors might be, we can 
calculate the mass bias from our mismatch calculations. For 
the least accurate case we have considered here, % — 0-5, with 
a hybrid matching frequency of Mco m = 0.05, the worst mass 
bias is 0.54% for M = 22 M . The results for masses up to 



100 M are shown in Fig. 1 1 We note for comparison that the 
results in Ref. fl23l suggest that the PN waveforms that can be 
used for searches for binaries with M < 1OM have relatively 
high mass biases, around 20% for M = 1OM . In addition, 
the results in Ref. Il49l suggest that the mass bias from using 
35PN TaylorTl inspiral templates to detect phenomenologi- 
cal waveforms is around 10%. These results of course also 
suggest that hybrids that are accurate enough for parameter 
estimation are also needed for binaries below 1OM . 

We emphasize that the true bias in the mass, when opti- 
mizing over all of the intrinsic and extrinsic parameters, may 
differ from that indicated in Fig. 1 1 While these results sug- 
gest that mass estimation errors will be very small with these 
waveforms, we defer a detailed analysis of parameter estima- 
tion errors to future work. 



IX. PARAMETER ESTIMATION 

Our focus so far has been on GW detection, and a full 
parameter-estimation study is beyond the scope of this paper; 
indeed such a study would require a complete waveform fam- 
ily, and would be more appropriately performed using phe- 
nomenological or EOB models. However, we can make some 
observations about parameter estimation. 



X. CONCLUSIONS 

In constructing hybrid PN+NR black-hole-binary wave- 
forms there are three evident sources of error: the error in 
the numerical waveforms, the error in the PN waveforms, and 
the error introduced in the hybridization process. The error 
measure relevant for GW detection is the mismatch error, and 
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in the context of current GW searches this error is required to 
be less than 3%. We know from previous work [22 29] that 
the mismatch error in the numerical waveforms is within this 
requirement by several orders of magnitude. We have shown 
in this paper that the mismatch error due to the hybridization 
process is also very low, at ~ 0.03%. The mismatch error due 
to the PN contribution to the hybrids, however, can be much 
larger, and dominates the error budget in all hybrids that can 
be produced with current PN and NR results. 

Motivated by these observations, we have attempted to ad- 
dress the question of how many NR cycles (i.e., how few PN 
cycles) must be used in a hybrid in order for the hybrid to 
be sufficiently accurate for GW detection purposes. Our ap- 
is to construct 



proach, which we justify in detail in Sec. VI 
TaylorT4+NR hybrids, and to treat these as the "true" GW sig- 
nal (NRl), and then to calculate the mismatch between these 
and hybrids of TaylorTl+NRL- The mismatch is optimized 
with respect to the total mass of the binary, which we argue in 
Sec. |VII| is close to the mismatch resulting from a calculation 
of the fitting factor, which is the relevant quantity for GW de- 
tection. This allows us to make the most conclusive statements 
possible with currently available numerical simulations. 

In the equal-mass nonspinning case we find that very few 
NR cycles are necessary; five cycles (two orbits) are sufficient 
to be within the 3% mismatch requirement. Phase compar- 
isons between PN and NR results suggest that the relative 
accuracy between TaylorTl and TaylorT4 does not change 
significantly as the mass ratio is increased to q = 4 ll22l . so 
we might expect that the length requirements are roughly the 
same for larger mass ratios. On the contrary, we find that the 
length requirements increase with the mass ratio, and for q = 4 
15 cycles (seven orbits) are necessary. We conclude that (1) 
we cannot infer NR-waveform length requirements directly 
from PN-NR phase comparisons, and (2) at higher mass ra- 
tios q > 4 we will require yet longer NR waveforms, but more 
extensive studies will be required to determine how many. 

Due to the greater PN errors for spinning binaries, more cy- 
cles are generally needed than in the nonspinning case. Our 
results are summarized in Tab. [II] We have considered only 
equal-mass spinning cases, but we can conclude from these 



results that for unequal-mass spinning binaries, more NR cy- 
cles will be needed than in any of the cases we have consid- 
ered here. 

A conservative summary of our results would be that in 
the cases we have considered, simulations of 10 orbits before 
merger should be sufficient to produce hybrids that are accu- 
rate enough for GW detection purposes. The resulting hybrids 
will be indistinguishable for SNRs in Advanced detectors of 
less than 30 for masses above 20 — 3OM . The important 
caveat to our results is that they apply only to the "cheapest" 
cases. Nonspinning binaries with q > 4 and unequal-mass 
spinning binaries will require more NR cycles, and we have 
not considered any cases with precessing spins. 

These conclusions depend strongly on the current state of 
the art. The advent of more accurate PN results (for example 
higher order spin contributions) could reduce these length re- 
quirements, as might a more robust quantification of the errors 
in PN methods. We have also not considered EOB results, for 
which the errors may be far lower, with a corresponding drop 
in the NR-waveform length requirements. 

In addition to being relevant to the construction of hybrids 
for analytic waveform models, our results also have direct 
bearing on two current efforts in the NR, data-analysis and 
analytical modeling communities. 

In the second stage of the NINJA project [ 16 1 hybrid wave- 
forms are being injected into detector noise to test search 
pipelines. The hybrids are constructed with matching fre- 
quencies in the range we have considered, but in most cases 
the matching will be around Mco — 0.07. In general our re- 
sults suggest that these hybrids will have a mismatch error of 
greater than 3% for masses lower than 20 — 3OM . 

The complementary NR-AR project [16| represents a 
community-wide effort to produce a large number of NR 
waveforms to calibrate analytic models. The nominal require- 
ment for these waveforms is that they include 20 GW cycles 
(10 orbits) before merger. From the point of view of con- 
structing hybrid waveforms for GW detection with typical PN 
approximants, our results suggest that these waveforms will 
be long enough along the branches of the parameter space 
that we have considered, but possibly not for high-mass-ratio 
cases, some unequal-mass spinning cases, and perhaps also 
cases with precession. 

We have presented a general method to assess NR wave- 
form length requirements, but our study only begins to address 
this question. Further work is needed to develop an easy-to- 
implement and robust method to estimate these requirements 
for configurations that have not yet been simulated in NR 
codes; work in this area is already underway l56l . It would 
be useful to make estimates of waveform length requirements 
for the construction of analytic models, both phenomenologi- 
cal and EOB; to extend this study to possible future detectors, 
like LISA 1571 and the Einstein Telescope ll58l ; and, finally, 
much more work is required to understand the length require- 
ments for parameter estimation, and to balance the needs of 
GW astronomy with the computational cost of numerical sim- 
ulations. 
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